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Abstract We identify coronal holes using a histogram-based intensity threshold- 
ing technique and compare their properties to fast solar wind streams at three 
different points in the hcliosphere. The thresholding technique was tested on 
EUV and X-ray images obtained using instruments onboard STEREO, SOHO 
and Hinode. The full-disk images were transformed into Lambert equal-area 
projection maps and partitioned into a series of overlapping sub-images from 
which local histograms were extracted. The histograms were used to determine 
the threshold for the low intensity regions, which were then classified as coro- 
nal holes or filaments using magnetograms from the SOHO/MDI. For all three 
instruments, the local thresholding algorithm was found to successfully deter- 
mine coronal hole boundaries in a consistent manner. Coronal hole properties 
extracted using the segmentation algorithm were then compared with in-situ 
measurements of the solar wind at ~1 AU from ACE and STEREO. Our results 
indicate that flux tubes rooted in coronal holes expand super-radially within 
lAU and that larger (smaller) coronal holes result in longer (shorter) duration 
high-speed solar wind streams. 

Keywords: Sun: Coronal Holes; Automated detection; Sun: Magnetic Fields; 
Sun: Solar Wind, Disturbances 

1. Introduction 

Understanding the connection between solar activity and space weather is im- 
portant for accurately predicting geomagnetic disturbances that can affect radio 
communications, satellite and spacecraft operations, ground-based power sys- 
tems and astronauts. Coronal mass ejections (CMEs) and flares are sources 
of energetic particles which can impact the Earth, and hence are of primary 
importance to space weather forecasting. As sources of high speed solar wind, 
coronal holes (CHs) are an important factor in producing recurring magnetic 
disturbances at Earth on time scales of days to months. It is well known that due 
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to the longevity of CHs, they can drive important physical processes on a larger 
time-scale than transient events like CMEs. Over its lifetime, a coronal hole can 
input as much energy into the heliosphere as a CME (Kavanagh and Denton, 
2007). Characterising CH properties allows us to determine the arrival time of 
high speed solar wind streams at Earth and other positions in the heliosphere. 

CHs are regions of the solar atmosphere that contain 'open' magnetic field 
lines which extend to interplanetary space and give rise to high-speed solar wind 
streams (Altschuler, Trotter, and Orrall, 1972; Krieger, Timothy, and Roelof, 
1973; Hassler et ai, 1999; Schwadron and McComas, 2003; Antonucci et ai, 2004; 
Fujiki et ai, 2005; Temmer, Vrsnak, and Veronig, 2007). The electron density in 
CHs is 2-3 times lower than that of the quiet Sun (QS), while the temperature is 
around 10^-10^ K (Munro and Withbroe, 1972; Wiegelmann, Xia, and Marsch, 
2005; Wilhelm, 2006). As a result of their lower density and emission measure, 
CHs appear as dark areas in X-ray and extreme ultraviolet (EUV) wavelengths 
(Reeves and Parkinson, 1970; Vaiana et ai, 1976). Magnetograms have also 
shown CHs to be of one predominant polarity (Wiegelmann and Solanki, 2004). 

There have been a number of attempts to study CHs from a theoretical per- 
spective. One of the main theories is the interchange model (Wang, Hawlcy, and 
Sheeley, 1996; LioncUo et ai, 2006), which is based on the assumption that the 
dominant process in the coronal 'open' field evolution is the reconnection between 
'open' and closed magnetic field lines. This view was initiated by observations 
showing streamer evolution (inflows and outflows) which could be explained with 
the closed flux opening into the solar wind and 'open' flux closing back into the 
streamer (Sheeley and Wang, 2002). Interchange reconnection could also explain 
the rigidly rotating CHs that are resistant to differential rotation (Timothy, 
Krieger, and Vaiana, 1975). The counter-argument to interchange reconnection 
is that the continuous opening and closing of coronal field lines would result 
in the injection of magnetic flux into the solar wind and cause a bi-directional 
heat flux (Antiochos et ai, 2007). The other main theory is the quasi-steady 
corona-wind model in which the coronal magnetic field is extrapolated from 
the photosphere using magnetograms. The simplest version, the source surface 
model, provides us with a topology of closed and 'open' flux regions by assuming 
a current-free corona and a completely radial magnetic field at a fixed surface 
(Luhmann et ai, 2003; Schrijver and DeRosa, 2003). To obtain the evolution of 
the coronal field, a sequence of the source surface solutions are used. In recent 
years, the quasi-steady model has been improved by including the solution of 
the full MHD equations, which made the assumption of a fixed surface and a 
current-free corona unnecessary (Pisanko, 1997; Rousscv et ai, 2003; Odstrcil, 
2003). 

To study their physical properties, CHs must be detected in a consistent 
manner. Initial attempts to extract coronal holes involved hand-drawn maps 
based on the two-day average of He i (10830 A) images and magnetograms 
(Harvey and Recely, 2002) made with the National Solar Observatory / Kitt Peak 
Vacuum Telescope (NSO/KPVT). In recent years, a number of automated selec- 
tion methods have been developed: Henney and Harvey (2005) used NSO/KPVT 
He I (10830 A) images, magnetograms and the above-mentioned hand-drawn 
CH maps; Malanushcnko and Jones (2005) used He i imaging spectroscopy to 
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determine CH boundaries; while Tonia and Arge (2005) have combined space- 
based EUV images with ground based He i (10830 A) and Ha observations 
and magnetograms to identify coronal hole regions. Chapman and Bromage 
(2002) and Chapman (2007) constructed synoptic CH maps using measurements 
made with the Solar & Heliospheric Observatory Coronal Diagnostic Spectrom- 
eter (SOHO/CDS; Harrisonl997). Delouille, Barra, and Hochedez (2007) used a 
multichannel segmentation method on EUV images to differentiate CHs, QS and 
active regions (ARs). More recently, SchoU and Habbal (2008) applied an image 
segmentation technique based on a combination of region and edge-based meth- 
ods using 171 A (Fe ix/x), 195 A (Fe xii), 304 A (He ii) images and Michelson 
Doppler Imager (MDI; Scherrer et ai, 1995) magnetograms to determine CH 
boundaries. 

Although determining the intensity threshold is of fundamental importance 
to accurately determine CH boundaries (and thus their properties), many of the 
methods discussed rely on arbitrary choices to determine the CH threshold. The 
main goal of our study was to produce a fast and robust algorithm to automat- 
ically identify coronal holes at any time of the solar cycle using a single image 
(195 A or X-ray) at a time. We introduced a three perspective identification 
of CHs using the Solar Terrestrial Relations Observatory/Extreme Ultraviolet 
Imager (STEREO-SECCHI/EUVI; Wuelser et ai, 2004) and the SOHO Extreme 
Ultraviolet Imaging Telescope (FIT; Delaboudiniere et al., 1995). The properties 
of several CHs were then compared with solar wind speed measurements made 
by the STEREO Plasma and Suprathermal Ion Composition instrument (PLAS- 
TIC; Blush et al., 2005) and the Advanced Composition Explorer (ACE; Chiu 
et al., 1998) to study the connection between the areas of CHs and solar wind 
flows at a number of perspectives at 1 AU. 

2. Observations 

The segmentation algorithm was tested on 195 A images made with the STEREO-| 
SECCHI/EUVI and SOHO/EIT instruments and X-ray images made with the 
Hinode X-Ray Telescope (XRT; Golub et al, 2007; Kano et al, 2008). The 195 A 
observations show the inner solar corona with a difference in the distribution of 
intensities between coronal holes and the quiet Sun. In this study, full resolution 
1024 X 1024 (FIT) and 2048 x 2048 pixel (EUVI and XRT) images were used. 
EIT 195 A images are normally produced every 12 minutes, and EUVI images 
every 10 minutes. The 195 A filter has a broad temperature response, which 
peaks at approximately 1.5x10^ K. 1024 xl024 pixel images made by XRT 
were also used in this study. The wavelength range of the XRT instrument is 
6-60 A, the temperature response of the Ti-Poly filter (27.4 A) used in this 
study peaks at around 9x10^ K. MDI line-of-sight magnetograms were used to 
differentiate coronal holes from filaments. 

The threshold of one large CH visible with the four satellites around 12-15 
January 2008 was studied in detail, while the long-term threshold variation was 
shown using observations obtained over 30 days in January 2008. We have also 
investigated the different line-of-sight magnetic field (i?LOs) distributions for low 
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intensity regions (LIRs) observed on 15 June 2008 and 22 March 2007. Finally, 
the characteristics of 46 large CHs were studied for a three- year period starting 
in 2006. Only major CHs with distinct solar wind profiles were included in the 
analysis to compare their area with high-speed solar wind stream durations. 



3. Data Analysis 

3.1. Local Intensity Thresholding 

Standard SolarSoft procedures (Freeland and Handy, 1998) were used to reduce 
and calibrate the images, which were then transformed from full-disk images 
to Lambert equal-area projection (LEAP) maps. The latter procedure projects 
a spherical surface onto a cylindrical surface whilst preserving the area and 
producing parallel latitudinal (and longitudinal) lines within approximately -60 
and -1-60 degrees latitude and longitude (McAteer et ai, 2005). The coordinates 
of the coronal hole boundaries are only reliable within the above mentioned 
range due to ncar-limb projection effects. 

Intensity thresholding is a powerful method for identifying features in images. 
It is well known that the intensity histogram of an image gives a multimodal 
distribution where each frequency distribution corresponds to a feature in the 
image (Gonzalez and Woods, 2002; Gallagher et ai, 1998; Gallagher et ai, 
1999; Belkasim, Ghazal, and Basir, 2003). In the case of 195 A images, each 
major feature (LIRs, QS and ARs) corresponds to a distribution in the intensity 
histogram. When an image is dominated by quiet Sun, the corresponding in- 
tensity histogram shows a unimodal distribution (top panel of Figure 1). When 
a relatively large part of the image is covered by LIRs, the resulting intensity 
histogram shows a bimodal distribution (bottom panel of Figure 1), where the 
local minimum corresponds to the boundary threshold. 

The minimum between the LIR and QS distribution was enhanced using a 
partitioning operation. By partitioning the solar image to sub-images, we ob- 
tained local histograms with more defined minima which aided the determination 
of thresholds. The sub-images used were 1/2, 1/3, 1/4 and 1/5 of the original 
LEAP map size. Each of the sub-image sizes were used to step across the map 
and obtain local intensity histograms from overlapping sub- images (Figure 2). 
This process provided more clearly discernible local minima for thresholding 
(Figure 3). This can be explained by different size sub- images having different 
proportions of QS and LIR and hence frequency distributions. The threshold 
minima are defined as local minima that are located within 30-70% of the QS 
intensity and are wider than 6 digital numbers (DNS). 

As mentioned above, the intensity range where the thresholds were identified 
depends on the QS intensity. For this reason the method is reliable at any time 
of the solar cycle regardless of the change in the overall intensity. It should be 
noted though, that with increasing solar activity, the acquired CH boundaries 
are less accurate due to bright coronal loops intercepting the line-of-sight and 
obscuring parts of the CH boundary. 
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Figure 1. Top: EUVI 195 A LEAP map showing no CHs, and the corresponding intensity 
histogram with a unimodal distribution. Bottom: LEAP map showing a largo CH, and the 
corresponding intensity histogram with a bimodal distribution. 
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Figure 2. The LEAP map is partitioned into overlapping sub-images. 



Once the final threshold was obtained, the LIR maps were created. Figure 4 
shows the same CH complex observed by four different instruments in two wave- 
bands (f95 A and X-ray). The observation times arc different to ensure that the 
CH is seen at the same central location observing with different satellites. We 
would like to draw attention to the large-scale similarity in CH morphology de- 
spite the use of different instruments and observation times. Also, the threshold 
histograms (right column of Figure 4) show how distinct the peak threshold is, 
which is then used to contour the LlRs. 
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3.2. Differentiating Coronal Holes from Filaments 

CHs were distinguished from filaments using MDI magnetograms. A magne- 
togram was found closest in time to the EUV observation and then differentially 
rotated to the EUV observation time. The LIR contours obtained from the EUV 
image were then overplotted on the differentially rotated magnetogram and the 
flux imbalance was determined for the contoured regions. CHs are well know 
to show one dominant polarity, while filaments have a bi-polar distribution of 
polarities (SchoU and Habbal, 2008). Figure 5 shows examples of LIRs observed 
on 15 June 2008 and 22 March 2007 respectively. The 15 June 2008 LIR is 
an example for a CH with a negative dominant polarity, while the 22 March 
2007 LIR is a CH with a positive dominant polarity (first and second row in 
Figure 5). The skewness of the field strength distribution was found to be -2.9 
and 6.4 respectively for the two CHs. The LIR highlighted on the bottom panel 
magnetogram in Figure 5 was found to be a filament based on the very small 
value of skewness (0.2) in its line of sight magnetic field strength distribution. 



4. Results 

Figure 6 shows a threshold consistency test. In this example, thresholds were 
acquired from 76 EIT 195 A images taken on 11 January 2008. It was expected 
that thresholds of similar intensity would be found due to the similarity of the 
images taken close in time (images were taken approximately every 12 minutes). 
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Figure 4. Left: STEREO A/EUVI, SOHO/EIT, Hinode/XRT and STEREO B/EUVI LEAP 
maps with LIR contours. The images show a complex CH structure approximately at the 
same location on the map. Right: Threshold histogram, where the peak value is the optimal 
threshold. 
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Figure 5. Left: EIT LEAP maps with LIR contours. Middle: MDl LEAP maps with LIR con- 
tours overlaid. Right: Bj^qs intensity histogram of the contoured regions. Regions highlighted 
in the first and second middle panels are CHs (Bj^qs histograms skewed), the third middle 
panel shows a filament {Bj^os histogram not skewed). 

The mean intensity threshold was found to be (17.6±1.2) DNS. Assuming that 
the intensity threshold distribution is normal, the true intensity threshold falls 
within ±14% of the calculated mean intensity threshold to 95% confidence. This 
confirms that the calculated thresholds are stable. Figure 7 shows the variation 
of the LIR threshold over an entire month (January 2008) using four instruments 
onboard STEREO, SOHO and Hinode. 

All coronal holes studied here were linked to significant rises in solar wind 
speeds at 1 AU. Figure 8 shows the CH positions at the start, maximum and end 
of the high-speed solar wind period on 14-24 June 2008. The images show that 
the west-most CH boundary was at approximately W50 when the high-speed 
solar wind stream arrived at ACE. The CH has passed over the west limb by 
the time the solar wind speed dropped back to normal levels (300-400 km s^^). 
The duration of the high-speed solar wind stream is related to the expansion 
factor of the magnetic flux-tube rooted in the CH. This is discussed in more 
detail below. 
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Figure 6. The threshold values obtained from 76 EIT 195 A images taken on 11 January 
2008. The gap in the data between 15:00 to 21:00 is due to a data drop-out. 

The relationship between this CH and its associated solar wind stream mea- 
sured at three distinct points in the heliosphere was explored using imaging 
and in-situ measurements from STEREO B, ACE and STEREO A. The top 
row of Figure 9 shows the CH complex on three days during 14-24 June 2008. 
The corresponding solar wind velocity measurements, given in the bottom three 
panels of Figure 9, show a high-speed solar wind stream sweeping around the 
inner hehosphere, from STEREO B to ACE to STEREO A. The high-speed 
stream was initially detected by STEREO B at 12 June 2008 19:10 UT, when the 
west-most boundary of the CH was at ~W53 from the perspective of STEREO 
B. Subsequently, the solar wind speed began to rise in ACE in-situ measurements 
two days later on 14 June 2008 11:00 UT. At this point in time, the CH's west- 
most boundary was at ~W50 as seen by SOHO. The high-speed stream finally 
traversed the point in the heliosphere monitored by STEREO A on 16 June 
2008 6:30 UT, at which time the western boundary of the CH was at '^W43. 
The complementary imaging and in-situ capabilities of SOHO, ACE and the two 
STEREO satellites offer a unique opportunity to study the relationship between 
CHs rotating across the solar disk and high-speed solar wind streams sweeping 
through the heliosphere. 

As indicated previously, it is well-known that there is a relationship between 
the size of CHs and the duration of their resulting high-speed solar wind streams, 
whereby larger CHs produce longer duration high-speed streams. We investi- 
gated this relationship using 46 CHs which were automatically identified using 
our segmentation algorithm. Figure 10 shows the resulting plot of high-speed 
solar wind stream duration as a function of CH area. The correlation coefficient 
obtained was ~0.6, which indicates a relationship between a CH area and the 
duration of a corresponding high-speed solar wind stream. 

5. Discussion and Conclusions 

We have developed a fast and robust CH identification technique that can be 
used with a single 195 A or X-ray image, in conjunction with a corresponding 
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Figure 7. Variation of the LIR threshold intensity over a month (January 2008) obtained 
using STEREO/EUVI B, SOHO/EIT, STEREO/EUVI A 195 A and Hinode/XRT X-ray 
images. 
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Figure 8. Top: EIT LEAP maps showing the CH position on the day of the start, the 
maximum and the end of the high-speed solar wind period. Bottom: The corresponding ACE 
solar wind speed measurements for the studied period. 



magnetogram. The segmentation algorithm first identifies LIRs using the local 
histograms of either 195 A or X-ray observations, and then exploits MDI magne- 
tograms to differentiate CHs from filaments by the flux imbalance determined for 
the identified regions. This was achieved assuming that CHs have a predominant 
polarity as opposed to filaments. Once the CH regions were determined, the 
position, area and mean B^os were calculated for multiple CHs selected in the 
period 2006-2008. 

The main limitation of the method results from near-limb effects which in- 
clude: line-of-sight obscuration in EUV and X-ray images; increase of noise in 
Blos measurements, and interpolation effects caused by the Lambert equal-area 
projection. These near-limb effects can cause incorrect boundary determination, 
the inability to determine a dominant polarity and incorrect area determination. 
Also, some CH boundaries located near the limb were found to have different 
threshold intensities which depended on the degree to which they were obscured. 
For this reason, only CHs within -60 and +60 degrees latitude and longitude were 
included in this analysis. 

From a scientific perspective, the main goal of this study was to link the on- 
disk properties of CHs to high-speed solar wind streams measured with ACE and 
STEREO near 1 AU. The start, maximum and end of a high-speed solar wind 
stream was determined from ACE and STEREO solar wind measurements and 
the corresponding CH locations were shown on EIT and EUVI 195 A images. 
Figures 8 and 10 demonstrated the obvious signs of non-radial CH flux tube 
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Figure 9. Top: EIT LEAP maps showing the CH position on three different days when the 
high-speed solar wind reached STEREO B, ACE and STEREO A. Bottom: The solar wind 
profiles obtained from the three satellites, with the line marking the start of the high-speed 
solar wind. 



expansion. It must be noted that the time for passage through a solar wind 
stream is hkely to be extended due to the dragging of the solar wind into the 
non-radial shape of the Parker spiral. We are currently extending this study 
to determine the expansion factor of CH flux tubes, and as a first step, have 
compared the area of 46 CHs to the corresponding high-speed solar wind periods. 
Although a large scatter was found in the data, a positive correlation was found 
(Figure 10). Initial estimates indicate expansion factors of approximately 2-10 
for flux tubes extending from the Sun to 1 AU. This assumes that the duration of 
a high-speed solar wind stream goes as Tduration ~ fi'^^-^CH^idth/SQQ), where / 
is the flux-tube expansion factor and CH^idth is the CH longitudinal width. This 
super-radial expansion of flux-tubes rooted in CHs is consistent with previous 
theoretical and observational results (Wang and Sheeley, 1990; Guhathakurta 
and Fisher, 1998; DeForest, Lamy, and Llebaria, 2001). 

In addition to CH flux tube expansion factors it is also important to de- 
termine when and where a particular portion of a solar wind stream emerged 
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Figure 10. Scatterplot showing the relationship between the area of CHs and the high-speed 
solar wind duration caused by the corresponding CHs. 46 CHs were included in this study. 



from the Sun. When attempting to determine when a high-speed solar wind 
stream will impact the Earth, the on-disk origin of the stream must be known to 
relatively high precision. Knowing the exact location of the CH when the high- 
speed plasma left, allows the determination of the path-length of the plasma. 
This provides the opportunity to study the modification of the Parker spiral 
through slow and high-speed solar wind interaction regions. Ultimately these 
methods will enable us to make more accurate estimates of the arrival times and 
durations of high-speed solar wind streams at the Earth and/or other points in 
the heliosphere. 

We conclude that an extensive study of the connection between the CH 
morphology, physical properties and the solar wind profile is needed to have 
a better understanding of coronal holes and their interplanetary effects. Our 
future plans are to extend our analysis for a solar cycle and link all the major 
CHs to solar wind speed rises. This will allow us to study the relationship of the 
CH properties and morphology with the corresponding high-speed solar wind 
profiles. 
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